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AUTOMATED RADIUS OF CURVATURE MEASUREMENTS 



RELATED APPLICATIONS 

Under 35 U.S.C. §1 19(e)(1), this application claims benefit of prior U.S. Provisional 
Applications No. 60/252,1 16, entitled "AUTOMATED RADIUS MEASUREMENTS" filed 
5 November 20, 2000 and incorporated herein by reference. 

TECHNICAL FIELD 

This invention relates to radius of curvature measurements, and more particularly to 
an apparatus and automated method for determining the radius of curvature for both concave 
and convex surfaces. 

5 10 BACKGROUND 

, : :j t 
.?;!.. \ 

A general method for determining the radius of curvature of an object's surface is 

|J based on a differential measuring technique. In the differential measuring technique, a 

63 

41 technician positions a reference element and the object to form a retroreflecting cavity by 

centering each element on an optical centerline and orienting the two elements so that they 

M 15 are parallel to each the other. A retroreflecting cavity forms when the two elements are 

spaced apart such that collimated light entering the cavity through the reference element is 

j*j brought to a focus on the surface of either the element or object. The null order of the 

U retroreflecting cavity characterizes how many times a light ray entering the cavity parallel to 

but off-axis with respect to the optical centerline, bouncing between the reference element 
20 and the object and intersecting the optical axis at either the reference element or object, 

strikes the surface of the object before it exits the cavity parallel to but off-axis with respect 
to the optical centerline. The null order changes as the object having a single radius of 
curvature is translated closer to or away from the reference element. The radius of curvature 
is calculated via a paraxial ray approximation by measuring the distance in cavity length 
25 between successive null orders. Differential techniques for measuring radii of curvature are 
discussed, for example, in "Differential technique for accurately measuring the radius of 
curvature of long radius concave optical surfaces" published in SPIE, 192, 75-84, 1979. 
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Typically, radii measurement techniques such as those described above required time 
consuming and error generating manual adjustments in determining the radius of curvature of 
an object. 



SUMMARY 

5 In general, in one aspect, the invention features an interferometric system for 

measuring the radius of curvature of a measurement object. The system includes a 
tunable coherent radiation source capable of emitting a radiant energy beam having a 
characteristic wavelength and of scanning the characteristic wavelength over a range of 
wavelengths; an unequal path interferometer which during operation includes a reference 
10 object and the measurement object and which receives a portion of the radiant energy 
beam from the tunable radiant energy source and generates an optical interference 
P pattern; a detecting system including a detector for receiving the optical interference 

yjj pattern; and a system controller connected to the tunable radiant energy source and the 

s'V 1 : 

gf detecting system. The controller is programmed to cause the tunable radiant energy 

63 1 5 source to scan the characteristic wavelength over the range of wavelengths while 

|2 concurrently monitoring the optical interference pattern via the detecting system and 

«i further programmed to calculate the radius of curvature of a surface of the measurement 

|t=:|t 

Li, object from the monitored optical interference pattern. 

Jj Preferred embodiments include one or more of the following features. The 

p 20 system controller is further programmed to calculate a phase shift in the measured optical 
interference pattern caused by scanning the characteristic wavelength. The system 
controller uses the calculated phase shift to compute the radius of curvature. The system 
controller is also further programmed to compute an optical path difference from the 
calculated phase shift and it uses the computed optical path difference to compute the 
25 radius of curvature. To compute the radius of curvature, the system controller is 
programmed to use numerical simulation or use recursive computation techniques. 
During operation the reference and measurement objects form ari interferometric cavity 
with the measurement object approximately in a null position and the system controller is 
further programmed to correct the computed radius of curvature to take into account 
30 possible initial positioning errors of the measurement object relative to the null position. 
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In some embodiments, the measurement object in the interferometric cavity is in a null 
position having a null order of n, wherein n is an integer at least as great as 2, and the 
system controller is programmed to compute the null order for the interferometric cavity. 

Also, preferred embodiments might include the following additional features. 
The tunable radiant energy source is a tunable light source. The unequal path 
interferometer is a Fizeau interferometer. The system also includes a wavelength change 
monitor configured to monitor changes in the characteristic wavelength of the radiant 
energy source. In that case, the system controller is arranged to receive output from the 
wavelength change monitor and is programmed to calculate the radius of curvature of the 
surface of the measurement object from a monitored change in the wavelength and the 
monitored optical interference pattern. In addition, the system controller is programmed 
to calculate the change in the characteristic wavelength and a phase variation of the 
optical interference caused by scanning the characteristic wavelength over the range of 
wavelengths and it is programmed to use the calculated wavelength change and the 
calculated phase variation to calculate an optical path difference in the interferometer, 
and then use the computed optical path difference to compute the radius of curvature. 

Also, the wavelength change monitor includes a detector and an unequal fixed 
path length interferometer (e.g. an optical fiber interferometer, a Michelson 
interferometer, or a Mach-Zehnder interferometer). The reference object has a curved 
surface (e.g. one which is concave). 

In general, in another aspect, the invention features an interferometric system 
includes a tunable coherent radiation source capable of emitting a radiant energy beam 
having a characteristic wavelength and of scanning the characteristic wavelength over a 
range of wavelengths; a wavelength change monitor configured to monitor a change in the 
characteristic wavelength of the radiant energy beam; an interferometer which during 
operation includes a reference object and the measurement object to form an interferometric 
cavity and which during operation also receives the radiant energy beam from the tunable 
light source and generates an optical interference pattern therefrom; and a detecting system 
comprising a detector for receiving the optical interference pattern. 

In general, in yet another aspect, the invention feature? a method for 
interferometrically determining a radius of curvature of a curved surface of a 
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measurement object. The method includes positioning the measurement object relative to 
a reference object to form an interferometric cavity; scanning a characteristic wavelength 
of a radiant energy beam over a range of wavelengths; recording an optical interference 
signal caused by the interferometric cavity as the characteristic wavelength of the radiant 
energy beam scans over a range of wavelengths; extracting a phase variation from the 
optical interference signal; determining a relative change in the characteristic wavelength 
which caused the extracted phase variation; determining an optical path difference 
between the reference and measurement objects from the phase variation and the relative 
change in frequency of the characteristic wavelength; and calculating the radius of 
curvature for the curved surface of the measurement object. 

Preferred embodiments of the invention might include one or more of the 
following features. The radius of curvature is calculated from the optical path difference 
and a null order of the interferometric cavity. The method also includes determining a 
null order of the interferometric cavity, e.g. by performing mechanical phase shifting. 
The operation of positioning of the measurement object involves positioning it relative to 
the reference object so as to form a retroreflecting cavity having a null order equal to n ? 
wherein n is an integer that is at least as large as 2. The method also includes extracting a 
2-dimensional phase variation from the optical interference signal and then using the 2- 
dimensional phase information to correct the calculated optical path difference between 
the reference and measurement objects. The operation of determining the relative change 
in the optical frequency of the characteristic wavelength involves monitoring the relative 
optical frequency change of the characteristic wavelength. 

The invention enables one to measure both concave and convex radii of objects 
spanning four orders of magnitude, e.g., about 10mm to about 150m, to high precision, e.g., 
0.001%. The invention also enables one to simultaneously provide a measurement of the 
surface profile of the object. 

Embodiments of the invention may include one or more of the following advantages. 
They enable an operator to determine the radii of objects automatically and without 
physically contacting the objects, i.e., non-contact. The automated operation of some 
embodiments of the invention facilitate measuring a high number of object radii such as in 
production environments. The invention enables one to build an instrument that is robust and 
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automatically corrects for variations in the initial positioning of the object within the 
instrument. It also enables one to rapidly determine the object radii, e.g. in less than one 
second. The invention also enables one to construct an instrument that is insensitive to 
temperature and Abbe errors. 

Other advantages of various embodiments of the invention include the ability to - 
measure the radius of curvature of spherical surfaces to high precision. Embodiments can be 
adapted to measure both concave and convex surfaces spanning a large range of curvatures, 
to require a minimum of manual preparation, to provide a means for automated measurement 
of similar surfaces in a production environment, to increase the throughput of such 
measurements, to concurrently measure test surfaces topography, and to minimize the 
influence of environmental effects that may degrade performance. 

The details of one or more embodiments of the invention are set forth in the 
accompanying drawings and the description below. Other features, objects, and advantages 
of the invention will be apparent from the description and drawings, and from the claims. 

DESCRIPTION OF DRAWINGS 

Fig. 1 is a schematic diagram of an apparatus for measuring the radius of curvature; 
Fig. 2 is an expanded schematic diagram of the apparatus of Fig. 1; 
Fig. 3 is a schematic diagram of a Fizeau interferometric cavity; 
Fig. 4 is an expanded schematic diagram of the Fizeau interferometric cavity of the 
apparatus of Fig. 2; 

Fig. 5 is a schematic view of the interferometer monitor of the apparatus of Fig. 2; 
Fig. 6 is a schematic view of the wavelength monitor of the apparatus of Fig. 2; 
Fig. 7 is a schematic view of a retroreflecting cavity for measuring convex surfaces; 
Fig. 8 is a plot of calculated contrast for coating reflectivities as a function of null 

order; 

1 dR 



Fig. 9 is a plot of the absolute value of sensitivity, 
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, as a function of a , F/R; 



Fig. 10 is a schematic diagram of an alternative wavelength monitor; 

Fig. 1 1 is a schematic diagram of a cavity geometry of Fizeau interferometer; 

Fig. 12 is another schematic diagram of a cavity geometry of Fizeau interferometer; 
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Fig. 13 is another schematic diagram of a cavity geometry of Fizeau interferometer; 
Fig. 14A is a plot of a intensity history distorted by multiple beam interference; and 
Fig. 14B is a plot of an undistorted intensity history. 
Like reference symbols in the various drawings indicate like elements. 

5 DETAILED DESCRIPTION 

An automated instrument provides a method for measuring both concave and convex 
radii of an object spanning over at least 4 orders of magnitude, e.g., about 10mm to about 
150m, to high precision, e.g., 0.001%. 

In general, the automated instrument includes a tunable light source, an unequal path 
10 interferometer cavity, and a detection system. The interferometric cavity is formed between 
a reference object having a known focal length and a curved measurement object. The 
Q curved measurement object is placed relative to the reference object in the interferometric 

% cavity to form a retroreflecting or a confocal cavity. In a retroreflecting or confocal cavity, 

s li:;f 

© an optical wave front returning from the cavity after one or more reflections from the 

15 measurement object surface, is substantially identical to an optical wave front reflecting from 

fj a surface of the reference object. 

s In the retroreflecting configuration, the beam comes to focus on a surface of either the 

hi* 

i !i it measurement object or the reference object. This is also sometimes referred to as the cats- 

H; eye configuration. The cavity has a null order that characterizes how many times a light ray 

p 20 entering the cavity reflects off the test object before it exits the cavity. The cavity null order 
changes as the measurement object is translated closer to or away from the reference object. 
Also, the cavity is characterized by a null position, which refers to the separation distance 
between the reference and test objects of the cavity with a particular null order. 

In a confocal cavity or confocal configuration, the wavefront is perpendicular to the 
25 surfaces of both the reference object and the test object. ... 

The cavity for a paraxial region produces an optical interference pattern having a 
single interference fringe when the tunable light source supplies coherent radiation into the 
cavity. At a fixed null position, the fringe intensity oscillates in a sinusoidal pattern as the 
tunable source is scanned over a range of wavelengths. The detection system monitors the 
30 oscillating interference pattern. The automated instrument measures the radius of curvature 
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of the measurement object, positioned relative to a reference object having a known focal 
length at a specific null position, by determining the optical path difference (OPD) of that 
cavity, i.e., the spacing between the reference object and the measurement object. As will be 
described in greater detail below, the instrument determines the OPD by scanning the tunable 
light source over a range of wavelengths, monitoring an oscillating interference pattern 
caused by scanning the wavelength, determining a phase variation of the monitored 
interference pattern as the wavelength is scanned, and determining the relative change in 
wavelength of light emitted by the tunable source that generated the phase variation. 

Referring to FIG. 1, a radius of curvature measurement apparatus 100 includes a 
tunable light source 10 that emits radiant energy 12 for an interferometer 30. In the 
described embodiment, the "radiant energy" is light but the term is meant to include 
electromagnetic energy of in many other frequency ranges. 

A wavelength monitor 20 measures the change in wavelength of light emitted by 
source 10 and an interferometer monitor 40 monitors the resulting interference produced by 
the interferometer. A central processing unit 50 is programmed to calculate the radius of 
curvature of a concave or convex surface of a test object 44 placed in interferometer 30 (See 
Fig. 2) based upon the null order of the interferometer cavity and on the optical path 
difference (OPD) of the cavity as derived from the electronic output of wavelength monitor 
20 and interferometer monitor 40. 

Large Radii (Radii of Curvature Greater than about 100mm) 

Referring to Fig. 2, tunable light source 10, e.g. a VORTEX™ series external cavity 
diode laser available from New Focus, produces a nearly collimated and monochromatic 
beam 12 of coherent optical energy that propagates through a beamsplitter 16. Beamsplitter 
16 redirects a portion of beam 12 into wavelength monitor 20, which may be a simple 
Michelson interferometer (Fig. 6), and allows the rest of beam 12 to propagate along optical 
path 14. An attenuator 1 8 positioned in optical path 1 4 after beamsplitter 1 6 controls the 
intensity of beam 12 propagating into interferometer 30. Interferometer 30 includes beam 
wavefront optics 22 and a Fizeau-type interferometer cavity 25. Beam wavefront optics 22 
include a diverting lens 26, a beamsplitter 28, and a collimating lens 32. Diverging lens 26 
converts beam 12 into a diverging wavefront 24 that passes through beamsplitter 28 and then 
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through collimating lens 32, which collimates diverging wavefront 24 to produce a piano 
wavefront 33. For ease of illustration, only a single ray of wavefront 33 is shown. Piano 
wavefront 33 propagates into interferometer cavity 25 formed between a surface 34 of a 
reference object 36 and a curved surface 42 of test object 44. The distance, L, by which 
5 reference object 36 and test object 44 are separated depends on the magnitude of the 
curvature of surface 42 and whether or not the surface is concave or convex. 

In the Fizeau-type interferometer cavity of Fig. 2, surface 34 of the reference object is 
flat and the curvature of surface 42 of the test object is concave. However, when measuring 
test objects having a convex surface and a large radius of curvature, the reference object will 
10 have a concave surface (See Fig. 3). The radius of curvature of the reference surface and the 
spacing between the two surfaces are selected to satisfy the relationship 0<F-L<R/2, 
where F is the focal length of the reference object, L is the gap between the reference and test 
O objects, and R is the radius of curvature of the test object. Mounts 38 and 47 hold reference 

object 36 and test object 44, respectively, and are translatable along optical axis 14 so as to 
W 15 change the distance between surface 34 and surface 42, i.e., change the OPD and null order 
g:;i of the cavity. As described below, apparatus 100 holds mount 47 stationary during radius of 

J curvature measurements whereas apparatus 100 may move mount 38, e.g., a Phase 

a Modulating Receptacle (PMR), to phase modulate the interferometric cavity and thereby 

14 determine the null order and generate a phase map of the curved surface of the test object, 

y 20 Referring to Fig. 4, as piano wavefront 33 impinges on reference object 36, a portion 

El of the wavefront reflects off of a partially reflective coating 34 ' on surface 34 as reference 

f " wavefront 48R. The appropriate reflectivity of coating 34' is determined below. The 

remaining portion of wavefront 33 propagates through coating 34' towards concave surface 
42 which, in turn, reflects a portion of wavefront 33 as a measurement wavefront 48M back 
25 towards reference object 36. Depending upon the radius of curvature of surface 42 and the 
OPD, a portion of measurement wavefront 48M will typically be reflected back and forth 
inside of the cavity at a specific number of locations on surface 42 and surface 34. At each 
point of reflection on surface 34, a portion of measurement wavefront 48M propagates 
through the reference object. Collimating lens 32 collects and focuses both reference 
30 wavefront 48R and measurement wavefront 48M into beamsplitter 28 which, in turn, reflects 
a portion each wavefront 48M and 48R towards interferometer monitor 40. 
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Referring to Fig. 5, interferometer monitor 40 includes a main aperture 52, a 
collimating lens 54, a beamsplittter 56, a charge-coupled device (CCD) camera 46 having a 
2-dimensional array of pixels, an aperture 58, a lens 62, and a cavity detector 64. Camera 46 
sends electronic signals to CPU 50 and a visual output monitor 49. Collimating lens 32 (Fig. 
4) focuses reference wavefront 48R and measurement wavefront 48M at main aperture 52 of 
interferometer monitor 40 so that only those portions of the reference and measurement 
wavefronts parallel to optical axis 14 are detectable by monitor 40. Collimating lens 54 
collimates the filtered components of wavefronts 48R and 48M onto beamsplitter 56 which 
transmits a portion of the collimated wavefronts into camera 46 and reflects a portion of the 
collimated wavefronts through aperture 58 and into lens 62. Lens 62, in turn, focuses the 
wavefronts onto cavity detector 64. Aperture 58 limits the solid angle viewing area so that 
detector 64 only receives the portion of the wavefronts near to the paraxial region. 

Referring to Fig. 6, wavelength monitor 20 includes an unequal path length 
interferometer 21 having mirror 66 positioned further from a beamsplitter 72 than a mirror 
68. A detector 74 records the interference image from interferometer 21. 

Referring back to Fig. 2, during operation, a system operator positions a concave test 
object in mount 47 of interferometer 30 to form a Fizeau-type interferometer retroreflecting 
cavity with the reference object. Observing the visual output of camera 46 on monitor 49, 
the operator translates mount 47 along optical axis 14 to find a null position of the 
retroreflecting cavity. A null position exists when the light beam intersects optical axis 14 at 
either surface of the cavity and the null order is defined as the number of times that the beam 
strikes the test object before being imaged through aperture 52 and onto the camera. The null 
order of the arrangement shown in Fig. 4 is two. 

The output of camera 46 includes a series of interference fringes which represents a 
2-dimensional interference signal between the reference and measurement wavefronts 
imaged into interferometer monitor 40. Camera 46 is a 12 bit high dynamic range camera 
which helps the operator to better visualize the fringe patterns from the interfering 
wavefronts. The number of interference fringes observable on the output of the camera 
depends upon whether or not the test object is positioned in a cavity null position. When the 
operator translates the test object into a cavity null position, the number of interference 
fringes is less than the number of interference fringes when the test object is not in a null 
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position. As a result, the operator observes the number of interference fringes on the camera 
while translating the test object to determine and locate the null positions. Once the operator 
nulls the cavity, apparatus 100 determines the radius of curvature automatically. 

Before initiating the automated measurement sequence, the operator can enter the 
specific null order, if known, into CPU 50. As will be described in more detail below, the 
CPU uses the specific null order and the measured OPD of the cavity in a recursion 
relationship to determine the radius of curvature of the test object. Using those same 
recursion relationships, the operator can also determine the specific null order by knowing 
both an approximate value for the radius of curvature of the test object and the distance 
between the test object and reference object. If an approximate value of the radius of 
curvature is unknown, apparatus 100 can automatically determine the null order by fixing the 
wavelength, translating mount 38 a known distance, and measuring the relative change in 
phase of the Fizeau interferometer cavity. 

Apparatus 100 measures the relative change in phase of the cavity via camera 46. 
Camera 46 sends electrical outputs to CPU 50 that are proportional to the intensity of the 
interference signal of the focused reference and measurement wavefronts on each pixel of the 
camera. As apparatus 100 changes the OPD of the cavity by moving mount 38, the 
interference signal from each pixel of camera 46 oscillates according to the relationship 

I o = C 0 + Ci cos(^) (1) 
where C 0 and Cj are constants, (p is the overall phase. A change in the overall wavefront 
phase is related to the change in the gap, L, between the reference and test object, and the 
null order via the expression 

A<p = 2nNkAL (2) 
where A<p is the change in overall phase; n is the index of refraction (typically equal to 1 for 
air); N is the null order of the cavity; k is the wavenumber equal to 2n/X, in which X is the 
center wavelength of beam 12; and AI is the relative displacement of mount 38. Using 
equation 2, CPU 50 calculates N by knowing how far it translated mount 38 and measuring 
the change in phase caused by translating the mount. . ~v - ^ 

CPU 50 also uses known-phase-shifting algorithms to analyze the oscillating intensity 
data for each pixel and thereby generate a 3-dimensional phase map of the test object's 
curved surface. CPU 50 stores the phase map into memory and, as described below, uses this 
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map to correct the calculated radius of curvature for errors which result from errors in 
positioning the test object at a cavity null. The phase map can also be used to generate a 
surface profile of the test object. The ability of apparatus to profile completely the test object 
depends on the area of the test object sampled by camera 46. For convex test surfaces, the 
sampling area increases as the null order increases. 

Recovering the phase from the oscillating intensity signal involves calculating the 
phase of the Fourier coefficient corresponding to the fundamental frequency of the intensity 
sine wave. This is done by Fourier transforming the signal, selecting out the complex 
coefficient a representing the amplitude of the frequency of interest and calculating 
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Phase shifting algorithms can be found, for example, in Optical Shop 



Testing (2 n Ed.), Chapter 14, p. 501 (Wiley, New York, 1992). Fourier-series phase shifting 
algorithms are discussed in Applied Optics 34(22), 4723-4730 (1995) by P. de Groot, and in 
U.S. Patent No. 5,473,434 and U.S. Serial No. 09/349,593 filed on July 9, 1999. 

Once the CPU obtains the null order either by input from the operator or by 
translating mount 38 and recording the resulting phase change, apparatus 100 holds mount 38 
in its original position and CPU 50 sends an electronic signal to light source 10 to scan the 
center wavelength of beam 12. As CPU scans the center wavelength, it concurrently records 
two series of intensity data (measurement points), one series measured by the interferometer 
monitor cavity detector and the other measured by the wavelength monitor, such that the data 
are synchronously related. This could be accomplished, for example, by assuring that each 
measured intensity measurement point resulting from a wavelength shift has a corresponding 
measured intensity resulting from a cavity phase shift. After completing the wavelength 
shift, the CPU analyzes the stored intensity histories from each monitor (monitors 20 and 40) 
with phase extracting algorithms, such as a sliding window FFT method, to determine the 
relative phase and wavelength shift between two widely separated measurement points in the 
shift sequence. Typically, the CPU determines the relative phase and wavelength shifts 
between an intensity measurement point recorded early in the shift sequence and an intensity 
measurement point recorded at a time near the end of the sequehceV Changes in the. center 
wavelength of beam 12 cause changes in the phase of the cavity according to 
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<P = (3) 

c 

where <p is the phase variation in the cavity, o is the frequency variation of beam 12, and D 
is the total cavity OPD. Apparatus 100 determines the optical frequency variation of beam 
12, 0 , by measuring the phase change in the wavelength shift monitor while concurrently 
monitoring the phase change in the Fizeau interferometer. The phase in the interferometric 
cavity of the wavelength monitor is 

<P m =k m L m (4) 

where <p m is the phase of the wavelength monitor cavity, k m is the wavevector of beam 12, 

and L m is the total OPD of the wavelength monitor cavity. Assuming the monitor cavity 

OPD is constant and using the expression for the wavevector, k=2nl\, k=v/c, and Eqn. 4, the 
optical frequency variation can be expressed as 

t> = " (5) 

2nL m 

Substituting Eqn. 5 into Eqn. 3 and rearranging provides 

4-L..D (6) 

9m 

D is approximately equal to 2NL, where L is the gap between the reference and test objects. 
Apparatus 100 calculates the OPD, D, with high precision by simultaneously measuring the 

(b c 

shift of the center wavelength, i.e., 0 - , and its resulting phase variation of the Fizeau 

interferometric cavity, i.e., q> = m 

c 

Once CPU 50 calculates the interferometer cavity OPD, it uses this value along with 
the null order to determine the radius of curvature of the test object. CPU 50 determines the 
radius of curvature using either a paraxial approximation or a numerical ray trace simulation. 
A paraxial approximation for concave surfaces is described by Gerchman and Hunter in SPIE 
Vol. 192 Interferometry p. 75-84 (1979) and results from applying a Gaussian image 
equation and a conjugate recursion formula. The paraxial approximation is further simplified 
by using Chebyshev polynomials which indicate that for concave mirrors, the gap between 
the reference optic and the test mirror, L , is related to the radius of test object, R , via 
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where N = 2,3... is the null order and the C N _ } are even order Chebyshev polynomials of 
the first kind having the form T N +i(x), (N= 1,3,5,...) where x = z N /R and of the second kind 
having the form Un(x), (N= 2, 4, 6, . . .) depending on whether N is even or odd. 

Referring to Fig. 7, the paraxial approximation for a convex surface having a radius R 
and a curved reference element having a focal length, F, is given by two relationships of the 
form: 

— = — + — , and — = — + — (8) 

r s, s; f r t t; 

where i = 0,1, denotes the bounce number off of either the reference element or the test 
object, the distances S i are measured relative to the mirror and the distances 7] are measured 
relative to the lens. For ease of illustrating the cavity, the measurement and test objects are 
shown as having planar surfaces. Using standard sign conventions, an equation relating S i 

and 7] reduces to an equation relating the radius of the mirror 7? to the gap L and the focal 
length of the lens F . For 4 th order (4 bounces off the test object), the system is symmetrical 
about the optical axis. The paraxial approximation for the upper half of the optical axis 
produces the equations 

2 112 11^211 , m 
— = — + — , — = — + — and — = — + — (9) 
R S 0 S' Q F T 0 Tl R S } S[ 

Following the sign convention (distances to the left of the reference element are negative and 
to the right are positive) results in four supplemental equations relating various distances: 

S;=-Z, r 0 =S;+Z, TZ = S } +L and S 0 =F-L (10) 
Substituting these relationships into the paraxial approximation to eliminate all the variables 
except L , F and R produces: 

A = ! + 1 (11) 

F R(F-L) [L RL iL 

2(F-L)-R 2L + R 
and solving for the correct root of R provides: - 

R ^ L{L- F)\%F - 1 61 + V32F] - (12) 

2F 2 +16L 2 -16LF 
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Table 1 lists the equations governing the dependence of the radius on the gap L (and for 
convex mirrors, also on the focal length of the transmission sphere F ) for the first 5 orders 
in the paraxial approximation. 
TABLE 1 



N 


R (Concave) 


R (Convex) 


2 


R = 2L 


R _-2L{F-L) 
F-2L 


3 


R=4L 


R _-4L(F-L) 
F-2L 


4 


R = 2(2 + 72)1 


_ -L(F-L)[4F-8L + 2yl2F) 
(F 2 -8FL + 8L 2 ) 


5 


R = 2(3 + y[5)L 


-L(F- L)(6F -161 + 2^f) 
(F 2 -12FL + 16L 2 ) 


6 


R = 4(2 + fi)L 


- 4L(F - L){2F -4L + >/3>) 
(F 2 -16FL + 16L 2 ) 



Alternatively, numerical simulations, such as ray tracing programs, can be used to 
calculate the radius of curvature. For instance, starting with a ray at the edge of the reference 
element, subsequent intersections of the ray with the test object and reference element can be 
calculated depending on the values assigned to the gap (L) between the reference element 
and the test object and the radius of curvature of the test object (R). By inverting this 
problem and using iterative function optimizing routines, the CPU determines the value of R 
that best fits a particular measured gap L. 

Once CPU determines a value for the radius of curvature, it corrects this value to 
account for errors resulting from the operator's positioning of the test object at a cavity null 
position. As discussed above, the operator uses the visual output from camera 46 to 
determine when the test object is approximately in a null position. In general, due to human 
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limitations, the operator cannot translate the test object to an exact null position and, as a 
result, positioning of the test object slightly behind or in front of the true null position causes 
fractional errors in the measured OPD and thereby produces an error in the computed radius 
on the order of about 0.03%. The positioning errors also produce focusing or imagining 
errors through the interferometer monitor which, in turn, result in an approximately quadratic 
phase deviation in the phase map of the test surface. Selberg discusses positioning and 
imaging errors in radius interferometry in Opt. Eng. 31, 1961-6 (1992). 

As described above, apparatus 100 uses camera 46 to record a 2-dimensional intensity 
history of the interfering wavefronts as CPU translates mount 38. The CPU compares the 
measured phase map to a theoretical phase map. The CPU generates the theoretical phase 
map using ray tracing routines, the measured OPD, known null order and focal length of the 
reference element (if applicable), and the calculated radius of curvature. The CPU iteratively 
changes the value of the measured OPD, calculates a new theoretical phase map, and 
compares that new map to the measured phase map. The CPU can compare the theoretical 
and measured phase maps by monitoring the difference in phase between the two maps as a 
function of 2-dimension position. The CPU repeats this procedure until it minimizes the 
difference between the measured and theoretical phase maps. Once minimized, the adjusted 
OPD value corresponds to an OPD value which would have been measured by apparatus 100 
if the test object were positioned to within about 1 micron of the exact cavity null position. 
CPU uses the null order and the adjusted OPD value, i.e., the value that produced the 
minimum difference between the measured and theoretical phase maps, to recalculate the 
radius of curvature. Apparatus 100 can also compare measured and theoretical 2- 
dimensional phase maps to correct errors in tip and tilt, as well as translation, of the test 
object. 

Once the system operator has positioned a first test object having a specific radius of 
curvature, apparatus 100 can be used in an assembly line to monitor the quality of optical 
components each manufactured to have the same radius of curvature. CPU 50 of apparatus 
100 can be programmed to translate mechanically successive optical components from an 
assembly line onto the preposition mount 47 and measure the radius of curvature of each 
element. „ : * -*v 
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Null Fringe Interference Visibility and Reference Object Coatings 

Null fringe visibility, or contrast, is important when the operator initially positions the 
test object in any retroreflecting or confocal position. In general, the reflectivity of coating 
34' on reference element 34 is selected such that null fringes are visible on monitor 49 over a 
large range of null positions by balancing the contrasts found in high and low order nulls. 
The visibility or contrast is given by 

VJ-^L. (13) 

It + Ir 

where Ir is the reference intensity and // is the test intensity. The reference intensity is 
defined as 

Ir = I 0 Rr 9 (14) 

where Rr is the reflectivity of coating 34', and the test intensity is defined as 

It = I 0 (\~Rr) 2 Rm N Rr N ' i . (15) 

where Rm is the reflectivity of the test object and N is the null order. Fig. 8 is a plot of 
calculated contrast for a coating reflectivities of 0.9, 0.8, 0.7, and 0.6 as a function of null 
order up to N=10. For N=2 to 10, a coating reflectivity of 0.7 provides good contrast both at 
low and high null orders. If the operator knows apriori where to position the test object 
relative to the reference object to obtain a specific null order, the operator can calculate 
which reflectivity provides the best fringe contrast and, based on this value, interchange 
reference objects into apparatus 100 having the desired reflectivity coating. 

Measurement Error 

Apparatus 100 measures three parameters, i.e., the null order, the wavelength shift 
and the cavity OPD. The null order is an integer and therefore is assumed to be determined 
by the CPU without any uncertainty. Uncertainty in the OPD measurements results from 
errors associated with monitoring the wavelength shift and the phase change of the 
interferometric cavity. Assuming uncorrected errors between the OPD and optical 
frequency (wavelength) shift measurements, the fractional error in the cavity OPD is: 

. . ; -J: : (16) 
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where a variable preceded by a delta sign represents the error in that variable. The fractional 
error in v is: 



Av 
v 





+ 








^ An > 



(17) 



Combining Eqns. 16 and 17 provides the total fractional uncertainty in OPD as 



AD _ 
^1 



\ L m j 



(18) 



The total fractional uncertainty in OPD, in turn, is used to determine the uncertainties 
associates with calculating the radius of curvature. 

For concave mirrors, the dependence of R on L as shown in Table 1 is linear. Using 
dR 

the relationship AR = — AL and since D = 2NL to an excellent approximation, AR can be 
dL 

rewritten asM = ■ For ^ < 10 the sensitivity, — — , is approximately 1 (as can be 

IN dL 2N8L 

verified via Table 1). As a result, for null orders less than 10, the uncertainty in the radius of 

curvature is approximately equal to the uncertainty in the OPD, i.e., AR « AD . 

AD dR 



For convex mirrors AT? , as in the concave case, can be rewritten as Ml xt = 

IN 



2NdL 



with the sensitivity 



dR 
8L 



obtained by using the equations in Table 1 . In general, the 



N 



sensitivity depends nonlinearly on both L and F . As an example, consider the case when 
N = 2. From Table 1, the sensitivity is: 

_ -2(F 2 -2FL + 2L 2 ) 
N=2 (2L-F) 2 

But at a cavity null, the gap is given by: 
F-R + ^R 2 +F 2 . 



dR 
dL 



(19) 



2 



(20) 



Substituting of (1 1) into (12) produces the following equation describing the sensitivity in 
terms of the ratio a = F/R ; 
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dR _ -2(q 2 +l-Vl + q 2 ' 

(i-VT^ 7 ) 2 

i a;? 



(21) 

#=2 



Referring to Fig. 9, a plot shows the absolute value of 



as a function of a for the 

N 



2NdL 

first 5 null orders. All sensitivities approach a finite limit for large a (F» R), but increase 

rapidly as l/a 2 for small a (F «R). At higher null orders, the sensitivity at F/R>1 also 

decreases from about 1 to about 0.1 with increasing null order. As a result, for F/R >1, the 

uncertainty in the radius of curvature is approximately equal to or less than the uncertainty in 

the OPD, i.e., AR < AD . 

The fractional error in the radius, AR/R, expressed as a percent of the true radius can 

be determined by the CPU and is useful in assessing the performance of apparatus 100. For 

concave objects, AR/R should decrease as R increases since AR « AD is essentially 

constant. For convex mirrors, the CPU calculates the fractional uncertainty in the radius by 

f . , . t ... I dR 
multiplying the sensitivity 



IN dL 



and a ratio AD/R. For F/R >1, the fractional 

N 

uncertainty decreases (increases) as R increases (decreases) because the sensitivity is 
essentially constant. For F/R<1, the fractional uncertainty increases as R increases because 
for fixed F , the sensitivity increases more rapidly than R . As a result, the fractional 
uncertainty in the radius for convex surfaces is minimized when F/R~\. 

Errors due to high order optical aberrations and imperfections in the optical 
components of apparatus 100 are expected to be smaller than the errors discussed above, but 
can be determined numerically. Departures from pure sphericity in the mirror surface also 
can be determined by the apparatus via an analysis of the intensity data detected by camera 
46. 

The wavelength monitor OPD stability effects the accuracy at which apparatus 100 
calculates the interferometer OPD. In general, longer wavelength monitor OPDs provide a 
more accurate measure of the wavelength shift but are more difficult to keep stable. A 
typical wavelength monitor has an OPD of at least about 250mm, is stable to better than 1 
part in 10 6 over an expected operating temperature range, and is insensitive to moderate 
acoustical and seismic influences. Wavelength monitor cavities exhibiting the desired 
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stability can be produced by assembling the interferometer on a bed of invar (r c = 8 -10" 7 ). 
If the temperature of the bed is known to 1°C then the fractional uncertainty in the cavity 
monitor length is « 10* 6 . If higher degrees of stability are necessary, the fractional 
uncertainty can be reduced by another factor of ten by replacing the bed of invar with 
Zerodur. Alternatively, a higher degree of stability can be achieved by using a fiber 
interferometer. 

Referring to Fig. 10, a fiber interferometer 200 includes a long length of optical fiber 
202 spiraled and embedded into a thermo-electrically cooled base 204. Two couplers 206 
and 208 separate and then recombine the fibers to produce long and short optical paths. The 
OPD of interferometer 200 is large, e.g., about 1m. Interferometer 200 is housed in a small 
package (approx 50mm X 50mm X 10mm). The thermal expansion coefficient of quartz, the 
material of optical fiber 202 and couplers 206 and 208, is less than 10~ 6 and the index 
temperature sensitivity is about 10" 5 . Assuming an effective wavelength monitor OPD 
temperature variation coefficient of fiber 202 of about 10" 5 and a temperature control of base 
204 to about 0.0 1°C, the cavity length can be stabilized to about 10" 7 . 

The detectors used in both the wavelength shift monitor and test cavity monitor are 
standard silicon devices. The necessary detector/amplifier bandwidth is determined by the 
maximum fringe frequency, / , which depends on the optical frequency sweep rate v and 
the longest anticipated test cavity OPD, Z) max , via 

f = ^L, (22) 

c 

For a cavity OPD of 20m, a 50GHz sweep and a sweep frequency of 100Hz, the maximum 
fringe frequency is only 333KHz. In general, the detector bandwidth must be at least 3 times 
this frequency, e.g., 1MHz, or higher. At a 1MHz sample rate, less than 20K samples of 
16bit data are acquired by each detector and sent to the CPU during the sweep. Both the 
wavelength monitor and the Fizeau interferometer detectors are sampled simultaneously 
thereby eliminating any requirements regarding how linearly the CPLf shifts the center 
wavelength of the light beam. As discussed above, the CPU determines the phase evolution 
in each cavity, i.e., wavelength monitor and Fizeau interferometer, from these data using a 
phase extraction algorithm such as the sliding window FFT method. 
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Small Radii (Radii Less than about 100mm) 

Using the cavity geometry of Fig. 2 for measuring concave surfaces with small 
radii results in fringe patterns on camera 46 that are too small to be reliably used for 
positioning the test object at a null position of the cavity. For convex surfaces, as the 
radius decreases below about 100mm, the amount of the surface area of the test object 
used in the radius measurement also decreases thereby reducing the area over which the 
apparatus can generate a phase map of the test object's surface. These problems can be 
addressed by using different cavity geometries as described below. 

As shown in Fig. 1 1, a cavity geometry 300 for measuring concave surfaces with 
small radii includes a concave reference element 330 and a test object 350 having a 
concave surface 355. Concave surface 355 is positioned relative to element 330 such that 
the distance, M, from the object's surface to focal point, P, of element 330 is equal to 
radius of object 350. Cavity geometry 300 produces a single null fringe on the CCD 
camera (not shown). The size of this fringe depends only on the relative numerical 
apertures (NA) of the reference element and the test object. A single large null fringe 
representing a wavefront produced from the entire surface of the test object is observed 
on the visual output of camera 46 provided that the NA of the reference element is larger 
than the NA of the test object. Using the visual output of the camera, the system operator 
translates test object 350 until a single large fringe appears which indicates that object 
350 is positioned relative to element 330 to create cavity geometry 300. 

As shown in Fig. 12, a cavity geometry 400 for measuring convex radii less than 
100mm has a null order of 1 . Cavity geometry 400 includes a concave reference element 
430 having a focus and a test object 450 having a convex surface 455 and a back surface 
460. The operator positions back surface 460 at the focal point, P, of element 430. Once 
positioned, apparatus 100 determines the radius of curvature as described above, except 
that the CPU does not need to determine the null order of the^cavity. The CPU translates 
mount 38, records intensity histories on the camera, calculates a measured phase map of 
the test object's surface, and stores the map in memory. Apparatus 100 measures the 
OPD by shifting the center wavelength of the light beam while simultaneously 
monitoring the change in phase in the wavelength monitor and Fizeau interferometer 
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monitor. The CPU uses these values to calculate a measured OPD and uses this value 
with the known focal length of the reference element to determine the radius of curvature 
of the test object. The CPU also corrects positioning errors by comparing a measured 
phase map to a theoretical map produced from the known focal length and measured 
OPD. 

Using the paraxial approximation, the measured OPD (L) will be equal to 2(F+R) 
(2(F-R)) for concave (convex) mirrors when correctly nulled. To a first approximation, 
the radii are given by (L-2F)/2 and ((2F-L)/2). The procedure for correcting the 
positioning errors presupposes that F is accurately known. In general, the focal length of 
a reference element is not constant but changes with temperature. Apparatus 100 
provides a convenient solution to this problem by calibrating the reference element focal 
length. 

Referring to Fig. 13, a system operator positions a curved surface 605 of a test 
object 610 at the focal point, P, of reference object 620, i.e., at the cats-eye position. 
Once in this position, the CPU shifts the center wavelength of the light beam and 
simultaneously records intensity data from the wavelength monitor and Fizeau 
interferometer monitor. Using this data, the CPU determines the OPD and calculates the 
focal length of the reference object. The CPU can also correct deviations from the 
perfect cats-eye position and the focal length by iteratively comparing the measured 
phase map at the cats-eye position with a theoretical phase map created with the known 
focal length and measured OPD. 

The radius fractional uncertainty, AR/R, depends linearly on the ratio F/R for both 
convex and concave small radius mirrors. Thus for fixed F the fractional uncertainty 
increases as R decreases. Assuming an OPD fractional uncertainty of 10' 6 , and including 
null cavity correction and phase errors, the fractional uncertainty for convex and concave 
surfaces linearly tracks the OPD error as 0.0001% per unit F/R. At F/R of 1, the 
fractional uncertainty is 0.0001% and at F/R of 10, the fractional uncertainty is about 
0.001%. 

Depending on the reflectivity of the reference and test surfaces, the confocal 
cavities shown in Figs. 1 1 and 12 act as high finesse Fizeau cavities in which there are 
multiple retroreflections which propagate into the interferometer monitor. The 
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reflectivity of the cavity R is related to the reflection amplitude for multiple beam 
interference by the relationship 



with p=kL, L is the OPD; and rl and r2 are the complex reflection amplitudes of the two 
surfaces in the cavity. For high surface reflectivities, the presence of multiple beam 
interference creates integer multiples of the fundamental frequency at which the intensity 
oscillates on the wavelength shift monitor and the cavity detector. The multiple beam 
interference distorts the intensity interferogram such that the intensity variation caused by 
moving the reference element or shifting the wavelength will no longer be a pure sinusoid 
(See Eqn.l). Fig. 14A shows the expected intensity variation for an aluminum mirror 
(80% reflectivity) and a reference surface made with a 70% reflective dielectric coating 
resulting from multiple beam interference. Fig. 14B shows the intensity variation for the 
retroreflecting cavity of Fig. 2 which does not suffer from multiple beam interference. 

In general, the presence of multiple beam interference is not a serious problem 
when determining the OPD. First, the quantity, quality and density of the intensity data 
recorded by the cavity detectors is high enough so that the band-limiting characteristics 
of the phase shifting algorithms used for extracting the phase in the OPD measurement 
easily track the fundamental frequency by reducing the sensitivity of the algorithm to 
higher-order frequencies caused by the multiple bearn intejference. Second, any phase 
error is cyclic and averages out over the multiple cycles recorded for the OPD 
measurement. ;; : - : : -/ r 7 . ... v. . V- - 

The multiple beam interference also creates potential problems in the intensity 
data recorded by the CCD camera for computing a phasg map and surface profile of the 
test object. In general, due to the limited operating rate of the CCD camera apparatus 100 
acquires significantly less data from the camera with sparse*- density than it does from the 
cavity detectors. As a result, the phase shifting algorithms cannot isolate the fundamental 
frequency as well as they did for the OPD measurement. For example, in a 13 frame 




(23) 



where 




(24) 
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algorithm, apparatus 100 records 13 intensity frames with 45 degree phase shifts between 
each frame. This particular algorithm has band-limiting characteristics which reduce the 
sensitivity of the algorithm to all even and odd harmonics of the fundamental frequency 
up to order seven but produce almost 20 nm of Peak-to-Valley (PV) error when used to 
5 analyze the intensity history shown in Fig. 14 A. One method for reducing this error 
involves reducing the magnitude of the multiple beam interference by changing the 
reflectivity of the reference object. For instance, if the reference object were left 
uncoated (a 4% reflectivity) the phase error associated with the 13 bucket algorithm 
drops to less than 0.2 nm PV. Of course, a lower reference reflectivity also reduces the 
10 contrast available on the CCD camera for aligning the test object and recording intensity 
data for phase mapping and surface profiling. The change in the contrast from a 4%-80% 
cavity to a 70%-80% cavity decreases by about a factor of 7, but still provides over 500 

p gray levels on a 12 bit CCD camera for the aligning and recording intensity data. 

^;;| Alternatively, apparatus 100 can shift the center wavelength of the light beam 

IB 15 rather than translating the reference object to acquire intensity histories on the CCD 

63 

camera. Shifting the wavelength generates a longer phase excursion, i.e., produces many 
more oscillations in the intensity history, than mechanically moving the reference object. 

ht'n 

r 

* Due to the increase in phase excursion, the CPU is able to record an increased number 

E j 

. it and density of intensity samples which, in turn, reduces the sensitivity of the phase 

H 20 shifting algorithms to non-fundamental frequencies. The CPU processes the intensity 
Q data by using a phase shifting algorithm, such as a windowed FFT, on each pixel's 

^ intensity history. For instance, apparatus 100 shifts the wavelength and the CCD camera 

records 64 consecutive frames at 22.5 degree phase increments. The acquisition would 
take only l A second for a 120Hz camera, and the distortion induced phase error for a 
25 70%-80% cavity would be less than 1 nm PV, i.e., improvement by a factor of 20. 

A number of embodiments of the invention have been described. Nevertheless, it 
will be understood that various modifications may be .qiade without departing from the 
spirit and scope of the invention. For example, although measuring small radii of 
concave surfaces are described in the context of cavity geometries having the concave 
30 surface at the first null position, small radii of concave surfaces can be measured via 
cavity geometries having hull positions greater than one. In these scenarios, the cavity 
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geometry requires a larger working surface, i.e., a larger OPD, relative to larger radii. In 
addition, the relative change in the characteristic wavelength emitted from the light 
source can be determined by monitoring electronic signals from the light source. For 
instance, the light source can be calibrated so that a specific change in a driver current 
corresponds to a known change in optical frequency output. Monitoring the relative 
change in driver current could, in turn, be used to determine the relative change in the 
optical frequency of the characteristic wavelength. 

Examples: Fractional Error Simulations 

The percent fractional error obtained from simulating measurements on concave 
and convex mirrors having radii ranging from 100mm to 150m are listed in Tables 2 and 
3. The wavelength shift detector was assumed to have 14 effective bits and cavity 
detector was assumed to have 12 effective bits. Results from a Sliding Window DFT 
simulation have shown that with this precision, the uncertainty in the phase evolution 

would be A<j> = 5X1 0" 5 and Acp = 2X1 0" 4 for a 0.3m wavelength monitor cavity and a 
lm OPD in a Fizeau interferometer cavity, respectively. Using these numbers and 
assuming a 50GHz optical frequency tuning range, the contribution to the fractional error 
in D from both phase measurements is less than 0.00002%. The fractional uncertainty in 
the optical length of the wavelength shift monitor was assumed to be 10* 6 (0.0001%) 
making it the dominant source in the fractional error in D. The results are displayed 
assuming null cavity error correction. Null cavity gap errors are assumed to be 1 micron 
after correction. 

The tables display all the features predicted by the paraxial analysis. For concave 
mirrors AR/R drops as R increases, and for convex mirrors, a minimum in AR/R is 
observed near F/R=l . In fact, if the null cavity "error can be ignored, the minimum 
occurs exactly at F/R=l for all orders. 
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Table 2: Simulated Errors for Measurements of Convex Radii 



Radius (m) 


Focal 
Length (m) 


Order 


F/R 


Gap (mm) 


- (%) 


0.1 


1 


5 


10 


910.49 


0.00113 


0.5 


1 


2 


2 


809.02 


0.00072 


1 


1.5 


3 


1.5 


911.44 


0.00025 


5 


4.5 


6 


0.9 


1002.88 


0.00007 


10 


4.5 


4 


0.45 


1011.99 


0.00021 


50 


4.5 


4 


0.09 


713.01 


0.00174 


150 


4.5 


4 


0.03 


676.25 


0.00568 



Table 3: Simulated Errors for Measurements of Concave Radii 



Radius (m) 


Order 


Gap (mm) 


R } 


0.1 


2 


49.94 


0.00199 


0.5 


2 


249.99 


0.00040 


1 


2 


499.99 


0.00020 


5 


4 


732.23 


0.00014 


10 


5 


954.91 


0.00011 


50 


11 


1012.68 


0.00010 


150 


19 


1022.90 


0.00009 



Accordingly, other embodiments are within the scope of the following claims. 
For example, though the wavelength monitor was described as a Michelson 
interferometer, other types of interferometers could be used, e.g. a fiber or bulk optic 
interferometer or a Mach-Zehnder interferometer. It should also be clear that the 
wavelength monitor only need s to be used if the wavelength uncertainty when scanning 
is too large. 
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